Amperean Pairing Instability in the U(l) Spin Liquid State with Fermi Surface and 

Application to « - (BEDT - TTF) 2 Cu 2 {CN) 3 
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Recent experiments on the organic compound K~(BEDT — TTF)2Cii2(CN)3 raise the possibility 
that the system may be described as a quantum spin liquid. Here we propose a pairing state caused 
by the 'Amperean' attractive interaction between spinons on a Fermi surface mediated by the U(l) 
gauge field. We show that this state can explain many of the observed low temperature phenomena 
and discuss testable consequences. 



The organic compound k-(BEDT-TTF)oCu2(CN) 3 
shows great promise as the first candidate 0, 01 which 
realizes the spin liquid state in dimension greater than 
one 3] . It is a quasi two dimensional material where 
each plane forms a half filled triangular lattice. While 
it is an insulator, there is no magnetic long range order 
and it behaves like a metal as far as the spin dynamics 
is concerned. The uniform spin susceptibility and the 
spin lattice relaxation rate \jT\T are finite in the zero 
temperature limit 0, 0|. Recently a specific heat mea- 
surement was reported which extrapolates to a linear T 
coefficient 70. The 7 and spin susceptibility forms a 
Wilson ratio close to unity. While an alternative inter- 
pretation in terms of Anderson localization has been pro- 
posed, we find that the variable range hopping fit to the 
resistivity when combined with the density of states 
derived from 7, imply a localization length of 0.9 lattice 
spacing. Such a short localization requires very strong 
disorder which make this interpretation implausible. 

Theoretically, the existence of the spin liquid state 
has been suggested from the studies of the extended t-J 
model|(| and the Hubbard model0,ll,l3 on the triangu- 
lar lattice on the insulating side of the Mott transition. In 
the spin liquid state, the low energy effective theory be- 
comes the U(l) gauge theory coupled with spinon which 
forms a Fermi surface 0. The spinon carries only spin 
but not charge and it contributes to the specific heat and 
thermal conductivity even in the insulating state. The 
system of the gapless spinon and the gauge field is a non- 
Fermi liquid statefnjllll an d exhibits singular tempera- 
ture dependence of specific heat coefficient, 7 ~ T~ 1//3 . 
However, recent specific heat measurement does not show 
this singular behavior ; 5j. Moreover, there exists a kink in 
the specific heat around 6K, which suggests that there is 
a peak in the electronic specific heat if the phonon part 
is assumed to be smooth. Around the same temperature 
the uniform susceptibility also shows a sharp drop be- 
fore it saturates to a finite value in the zero temperature 
limit. These suggest that there is a phase transition or 
a cross-over from the high temperature phase which is 



described by the non-Fermi liquid state to a low temper- 
ature Fermi liquid state. In the present paper we propose 
that the low temperature phase may be understood as a 
novel paired state of spinons that arises out of the U(l) 
spin liquid state with spinon Fermi surface. 

Inspired by Ampere's discovery that two wires carry- 
ing parallel currents attract each other we note that 
the interaction is attractive when the two spinons have 
parallel momenta. We therefore explore the possibility 
of pairing two spinons on the same side of the Fermi sur- 
face. Within a simple mean field treatment we demon- 
strate such a pairing instability of the spinon Fermi sur- 
face. The resulting state has a number of properties that 
are attractive for an explanation of the experiments in 
k - (BEDT - TTF) 2 Cu 2 (CN) 3 . In particular the pair- 
ing gaps out the gauge field so that the unpaired por- 
tion of the Fermi surface gives a linear specific heat at 
low temperature. We discuss various consequences of our 
proposal that may be tested in future experiments. 

Consider the system of spinon with Fermi surface in- 
teracting with (non-compact) U(l) gauge field in 2+1D, 
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C = ^(9o-i0-M)^ + ^C(-*V-a) 2 -0<T + ^/ l j7r 

(1) 

Here Xq is the imaginary time and x = (%i,X2) is the 
2d spatial coordinate. i/j a is the spinon field with spin 
a and /z, the chemical potential. Repeated spin indices 
are summed, oij = (4>, a) is the U(l) gauge field with 
« = 0,1,2, fij, the field strength tensor and g, the gauge 
coupling. We expect g 2 to be proportional to the charge 
gap and will ignore the last term in Eq. Q in the follow- 
ing. We choose the Coulomb gauge where V ■ a = 0. We 
are interested in the stability of local Fermi surface in the 
momentum space and we focus on a patch of Fermi sur- 
face which is centered at a momentum Q with |Q| = fop- 
Therefore we integrate out the spinon fields except for 
those in the patch. The massless spinons screen the tem- 
poral gauge field <j>. However, the transverse gauge field 
is not screened and it can mediate a long range inter- 
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action between spinons. The dressed propagator of the 
transverse gauge field is given by 
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where k = (fco,k) is energy- momentum vector and 7 = 
and Xd — \2-Kfn are ^he Landau damping and dia- 
magnetic susceptibility respectively, vp is the Fermi ve- 
locity and to, the mass which are averaged over the Fermi 
surface which has been integrated out. The transverse 
gauge field mediates an interaction between spinons 
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where V is the volume of the system, j3 = l/(fcsT) and 
to, the mass of the spinon in the vicinity of Q on the 
Fermi surface. It is noted that to is generally differ- 
ent from m if the Fermi surface is not perfectly spheri- 
cal. Motivated by the Amperean attraction, consider the 
pairing of two spinons with energy- momenta, p\ = Q+p, 
P2 = Q — p, where 2Q is the net energy-momentum of the 
pair with Q = 0, |Q| = kp and p, the relative energy- 
momentum with |p| << |Q|. Note that the pair is made 
of two spinons on the same side of the Fermi surface and 
it carries a large net momentum of 2k p. We decompose 
the two body interaction into pairing channel by intro- 
ducing the Hubbard Stratonovich field Ap " , 
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-r-D{k) and we used (Q + p) x k 



where v(k) 

Q x k. Pairing may occur in the singlet channel, i.e. 
Ap- 1 = A p , A p T = -A_ p and Ap T = A^ = 0, in which 
case A p is an even function of p, or the triplet channel 
where A p is odd in p. Now we integrate out the rest 
of the spinon field to obtain the Landau- Ginzburg free 
energy density, 



/[A] = A t (w-wnw)A + 0(A 4 ). 



(5) 



Here every product is a contraction of energy-momentum 
indices with a measure ^ . A is a vector with component 

A p , and v, n are matrices with elements v > = v(p —p) 
and U p > p = Vf3g(Q + p)g(Q - p)6 p > p . g(p) is the spinon 
propagator given by g(p) = -, — — - — ^ ; — ry- — with 

" J aV ' i{pa+MPa\ 2/3 sgn(p ))+e-p 

and e p , the spinon energy dispersion. 
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A ~ /J 2/3 1/3 

Here Vp is the Fermi velocity at the patch which generally 
differs from the averaged one (vf). 



The system is unstable against developing pairing am- 
plitude when an eigenvalue of the kernel (v — vILv) be- 
comes zero or negative. Defining $ = vA, we write the 
eigenvalue equation 



Epair® 



(6) 



where E pa i r is the eigenvalue and $ is the eigenvector. 
Along the direction of the eigenvector, the free energy 
density becomes /[A] w A^l-vII)* = (1 - E pair )A^vA 
and the system becomes unstable if E pa i r > 1. The 
components of 4> has the unit of energy and Eq. © 
is nothing but the linearized self-consistent equation for 
the anomalous self energy which gives rise to the gap in 
the quasiparticle spectrum. At zero temperature and in 
the thermodynamic limit, the matrix equation Eq. © 
becomes an integral equation, 



Epair® (p) 



dp 
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Hp - p )g(Q + p )g(Q - p )®(p )■ 

(7) 

First, we approximately solve the equation analytically 
by guessing an Ansatz for the eigenvector. Then we will 
check the validity of the analytic solution by solving the 
equation numerically without assuming a specific form of 
the eigenvector. We first consider singlet pairing. 

In order to guess the form of the eigenvector, we de- 
termine the important region of integration for p in Eq. 
(J7J. If the pairing interaction were instantaneous and 
the spinon did not have the frequency dependent self 
energy correction, the p Q integration would impose the 
constraint that both of the constituent spinons of a pair 
should be on the outside of the Fermi surface, that is, 

\ Vp P\\\ < fm' where p., (p ± ) is the momentum along 
(perpendicular to) the Q as is shown in the bottom inset 
of Fig. n In the presence of the frequency dependent 
interaction and spinon self energy, the sharp constraint 
is smeared out. However, the dominant contribution of 
the momentum integration still come from the region 

|i>fP||| < 2m which is denoted as the shaded area in the 
bottom inset of Fig. ^ This has been checked by per- 
forming a numerical integration of p . Knowing the im- 
portant region for p we consider an approximate Ansatz, 
<J>(po,pj_,P||) = ®(Po,P±)@(p±/m - \v F p\\\), where we 

2 

take the range of py twice larger than |i>fP||| < 2m m or ~ 
der to take into account the smearing effect. This ansatz 

is singular at the curve pp = which includes the point 
P|i = p± — 0. A better treatment will smear out this sin- 
gularity. For this Ansatz, the typical momentum transfer 

' k 2 

k = p — p also satisfies the condition |fc||| < << fcj_ 
and we can ignore the fen dependence in the gauge prop- 
agator. We can also replace ^ffi by vp because k is 
almost perpendicular to Q. We perform the fcy integra- 
tion in the Eq. JJJ and we obtain 
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x$(p + k ,p± + k±). 
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Here we use the simpler form of gauge propagator which 
is obtained from Eq. J5J in the limit w_F|k| >> |fco| and 
we keep only the leading frequency dependent term in 
the spinon propagator. Therefore, the integration for 
the energy /momentum should be understood as hav- 
ing a ultraviolet cut-off of the order of the Fermi en- 
ergy/momentum. 

It is noted that the right hand side of Eq. (JSJ 
is smooth as a function of po and depends on pa 
very weakly. Therefore we ignore the po dependence 
in the kernel and consider a frequency independent 
eigenvector. The gauge propagator 7o |fc |+x!|fcj :i is 

sharply peaked at fcj_ ~ (7o|fco|/Xd) 1 ^ 3 as a func- 
tion of k± and can be approximated by a delta func- 
tion (7oXdl fc o|)" 1/3 E s =±i^( fc ± - s(7oN/Xd) 1/3 ) and 
we can perform the k± integration. Changing the in- 

tegration variable k by t — s —kn we obtain the 
eigenvalue equation 
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dt 
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where A = [2m\(xd/jo) 2 ^ 3 ] is a dimensionless con- 
stant. If we consider the spinon pair right on the Fermi 
surface (p± = 0) and use the Ansatz $(pj_) = const., the 
right hand side of Eq. 10 is logarithmically divergent. 
This signifies that we can find an eigenvector which has 
an arbitrarily large eigenvalue. However, the momentum 
independent Ansatz cannot satisfy the eigenvalue equa- 
tion because the kernel strongly depends on p± . In view 
of the singular dependence of the kernel on p± , we con- 
sider an Ansatz 



\P±\ C 



(10) 



where a should be smaller than 3/2 in order for the 
eigenvector to be normalizable. This Ansatz solves 
the eigenvalue equation with the eigenvalue E pa i r ~ 

,„ 6 U In ( 1 + .? 4 > c 2i ) - , where c = measures the lo- 

(27z) z c V 3c^ + l J a ' mvp 

cal curvature of the Fermi surface. For small enough a 
the eigenvalue can be arbitrarily large. Thus within the 
present mean field treatment there is a pairing instability. 

Now we check the validity of the analytic solution by 
solving the eigenvalue equation numerically. We do not 
assume a specific form of $(p) in Eq. Q. Then the 
natural cut off for the fcy = p|| — p^ integration is k± 

not k\ I m because the coupling becomes small 



for fey > k±. Ignoring the k\\ dependence of the gauge 
propagator, we can cast the equation into a 2D integral 

equation by applying ^-g(Q + p)g(Q - p) on both 
sides of the Eq. Q. The resulting equation involves 
only two integrations and one can easily diagonalize the 
kernel M(p,p ) numerically to find the eigenvalue and 
the eigenvector. The largest eigenvalue corresponds to 
singlet pairing (even A p ) and as shown in Fig. ^ in- 
creases logarithmically with increasing L, where L deter- 
mines the mesh of the discrete energy and momentum as 
Afc = 2-k/L. The eigenvalue will become larger than 1 
for a large enough L and there exists pairing instability 
in the thermodynamic limit. The infrared divergence of 
the eigenvalue in the thermodynamic limit is consistent 
with the analytic result that the eigenvalue diverges as 
a — ► 0. Although not shown here, the numerically calcu- 
lated eigenvector is qualitatively consistent with the an- 
alytic Ansatz with a < 1. The second largest eigenvalue 
corresponds to triplet pairing and is also logarithmically 
divergent with a slope 10 times smaller than that shown 
in Fig. In the rest of the paper we assume singlet 
pairing, even though we should be mindful that triplet 
pairing is also unstable and may be preferred by short 
range repulsion. 

The origin of the mean field pairing instability should 
be contrasted with conventional superconductors where 
electrons with momenta p and -p form a pair which 
uses the whole Fermi surface to lower its energy. In the 
present case, spinon pairs carry a momentum of the or- 
der of 2kF- While the LOFF state also carries finite 
momentum 13], our case is fundamentally different be- 
cause pairing fermions on the same side of the Fermi 
surface severely restricts the available phase space. Con- 
sequently, the pairing instability found here can happen 
only if the pairing interaction is sufficiently singular. 

Now consider possible instabilities in the particle-hole 
channel. In the leading order of perturbation, the 2kp 
vertex function is logarithmically divergent if both the 
external momenta are on the Fermi surface and energy 
transfer is zero 14] . The singular vertex function en- 
hances the susceptibility of the spin density wave with the 
momentum 2kp. A mean field treatment similar to that 
above reveals the existence of 2kp density wave instabil- 
ities in both the singlet and triplet channels, albeit with 
non-trivial momentum dependence for the internal wave- 
function for the particle- hole pair. We emphasize that 
the mean field instability occurs in a channel where the 
eigenvector has a specific momentum dependence. The 
system can remain stable in other channel. We point out 
that the 2kp instability preserves the U(l) gauge struc- 
ture. Furthermore the Amperean pairing is favored for 
large local curvature c while the 2kp instability prefers 
small c. 

Theoretically the mean field results above should be re- 
garded as merely suggestive of possible low temperature 
instabilities of the spinon Fermi surface state. Lacking a 
better theoretical treatment we will take experiments as 
a guide for further discussion. 
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The NMR measurement does not observe the line 
broadening expected for the incommensurate spin den- 
sity waveplllfi|. Thus we discard the triplet spinon den- 
sity wave. The singlet spinon density wave state will have 
a non-Fermi liquid specific heat unless a further Am- 
perean pairing instability develops at low temperature 
to restore Fermi liquid behavior (see below). However in 
the latter case two separate transitions would have been 
expected as a function of temperature (for instance as 
visible signatures in the specific heat) which is not ob- 
served. Therefore, we focus on the scenario where only 
the spinon pairing occurs and explore some consequences. 

First, in the paired state gaps will open on the patches 
of the Fermi surface where the pairing occurs. The mo- 
mentum point at which the pairing instability first oc- 
curs depends on the details of the Fermi surface. In 
general there will be a number of preferred points re- 
lated by hexagonal symmetry. Once the pairing occurs 
on parts of the Fermi surface, the U(l) gauge group is 
reduced to Z 2 . Since the Z 2 gauge field is gapped, the 
low energy theory becomes the Fermi liquid theory and 
the remaining Fermi surface can remain gapless without 
further instability. This is consistent with the observa- 
tion that there exists a finite specific heat coefficient 7 in 
the zero temperature limit rather than the singular T~ 1//3 
behavior. The proposed spinon pair state will generically 
break lattice translation, rotational or even time reversal 
symmetries. As shown in the top inset of Fig. ^ sup- 
pose the pairing occurs at two distinct favored momenta 
Qi and Q2 with Ai, and A 2 the corresponding pairing 
order parameters. These are of course not gauge invari- 
ant but the gauge- invariant combination (Ai)*(A 2 ) is at 
non-zero momentum 2(Q 2 — Qi) and is also condensed 
as Ai and A 2 are individually condensed. This repre- 
sents a spontaneous breakdown of lattice symmetry. As 
we are discussing a spin system, this order corresponds 
to an incommensurate version of the valence bond solid 
(spin Peierls state) . However we emphasize that this bro- 
ken lattice symmetry state coexists with fractionalized 
spinons. The broken lattice symmetry implies a finite 
temperature phase transition. However due to the in- 
commensurate ordering the transition should be 2d X-Y 
like and shows no observable singularity in the specific 
heat. The translational symmetry breaking should cou- 
ple to lattice distortion and may be observable by X-ray 
scattering. 



A key prediction is that the low temperature thermal 
conductivity n ~ T like in a metal in contrast to the 
vanishing thermal conductivity expected in an Anderson 
insulator or the enhanced k ~ Ta for the spinon Fermi 
surface state with a gapless U(l) gauge field[8j. 

One may think that the sharp drop of the uniform 
spin susceptibility below 6-ftTllJ can be explained from 
the reduction of the spinon density of state caused by 
spinon pairing. However, this explanation is not cor- 
rect for the following reason. In contrast with BCS 
theory but in common with the LOFF state, the Am- 
perean pairing is not destroyed by the Zeeman limiting 
field because the spinon with up spin with momentum 
IQtI = IQI + HbH/vf and the spinon with down spin 
with momentum |QJ = |Q| — [IbH/vf can both be on 
the Fermi surface and paired without the energy cost of 
the Zeeman energy. This property is crucial in explain- 
ing the lack of field dependence up to 8T in the specific 
heat|5j- However, it follows that the uniform susceptibil- 
ity should not be affected by the opening of the pairing 
gap. We suggest that the reduction of the susceptibility 
may come from a contribution of the gauge field which 
is gapped out at low temperature. Details will be dis- 
cussed elsewhere fl6j . The spinon pairing state is not a 
superconductor because the spinon does not carry charge. 
However, if the charge gap is suppressed by driving the 
system across the Mott transition point with pressure 0, 
Bose condensation of the charge degrees of freedom con- 
verts the Amperean pairing state to a real superconduc- 
tor. One signature of this unconventional superconduc- 
tor is that the Knight shift will hardly change across 
the transition temperature, which is highly unusual for 
singl et p airing. This signature is consistent with recent 
data[l7|. 

In conclusion, aside from intrinsic theoretical inter- 
est, our proposal of a novel spin liquid state with paired 
spinons explains many of the unusual low temperature 
behaviors in k - (BEDT - TTF) 2 Cu 2 (CN) 3 and is 
amenable to further experimental tests. 
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FIGURE CAPTIONS 

Fig. 1 (color online) The largest eigenvalue of M(p,p ) as 
a function of the system size L where the mesh 
of the discrete energy and momentum is given by 
Ak — 2tt/L. Top inset : schematic picture of par- 
tial gapping of the spinon Fermi surface. Bottom 
inset : definition of p\\ and p±. 



